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' Abstract 

^ ; 

The momentum distributions of partons in bound nucleons are known 
Qh. to depend significantly on the size of the nucleus. The Fourier trans- 

^ I form of the momentum (xbj) distribution measures the overlap be- 

tween Fock components of the nucleon wave function which differ by 
a displacement of one parton along the light cone. The magnitude of 
^ I the overlap thus determines the average range of mobility of the par- 

' ton in the nucleon. By comparing the Fourier transforms of structure 

functions for several nuclei we study the dependence of quark mobil- 
ity on nuclear size. We find a surprisingly small nuclear dependence 
(< 2% for He, C and Ca) for displacements t = z ^ 2.5 fm, after which 
a nuclear suppression due to shadowing sets in. The nuclear effects 
observed in momentum space for xbj ~ 0.4 can be understood as a 
reflection of only the large distance shadowing in coordinate space. 



1. Introduction 

The difference between the quark structure functions of nuclei and those of 
free nucleons, first observed in 1982 by the EMC collaboration has gen- 
erated considerable experimental and theoretical interest p[. The measured 
nuclear structure function gives direct information about how the momen- 
tum distribution of quarks in nucleons is modified by nuclear binding effects. 
Several models have been proposed for such modifications, many of them 
based on a picture of the nuclear wave function in coordinate space, from 
which modifications of the momentum space distribution are surmised. 

The inclusive lepton scattering measurements determine the single parton 
distributions, but do not constrain parton-parton correlations in a model- 
independent way. Hence the inclusive data can only partially constrain the- 
oretical models for the nuclear effect. 

Here we wish to study the nuclear effects in coordinate space. The re- 
lation between quark distributions in momentum and coordinate space has 
been known for a long time, and involves no further model dependence than 
is needed for the usual QCD interpretation of the experimental structure 
functions ^, The phenomenological discussions of nuclear effects on 
parton distributions have nevertheless concentrated almost uniquely on mo- 
mentum space (see, however, Ref. While the momentum and coordinate 
space descriptions are in principle equivalent, insight into the physical mech- 
anisms may benefit from viewing the quark distributions in both spaces. For 
example, the standard explanation of the small x shadowing effect as due to 
the soft scattering on the nucleus of qq pairs created by the virtual photon 
well upstream of the target is most naturally discussed in coordinate space 



The cross section of deep inelastic lepton scattering (DIS) as a function 
of the photon virtuahty —Q^ and the Bjorken variable Xbj = x = Q"^ /2mv, 
where m is the nucleon mass and v is the energy of the photon in the target 
rest frame, can be parametrized as 

da _ Ana'^ F2{x,Q^) 

in terms of the structure function F2(x, Q"^) and the ratio -R(x, Q"^) = ail ax 
of the cross sections induced by longitudinally and transversely polarized 
virtual photons (in this formula, y = v/E is the fraction of initial lepton 
energy carried by the photon). Measurements on a variety of nuclear targets 
A have shown that the A-dependence of R{x,Q'^) is weak 0, so that the 
ratio of F2 structure functions for different targets Ai,A2 is approximately 
given by the ratio of measured cross sections: F^^/F,^^ = da^^ / da^"^ . The 
dependence of F^^ IF2^ is also known to be weak and will not be 
discussed below. 

According to perturbative QCD at lowest order in a^, the F2 structure 
function is given by 

F2(x, Q2) = ^ e\x\qlx, Q') + Q^)] (2) 

i 

where Q^) is a quark distribution in momentum space, i.e. the proba- 
bility that a quark of flavor i (having electric charge Cj in units of e) carries 
a light-cone fraction x of the nucleon momentum. 

The deep inelastic eN e'X cross section is related to the forward 
'j*N 'j*N hadronic matrix element 

T,. = J d'yeM^q ■ y){P\T[j,{yyUmP) (3) 
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where Fi = F2/2X at lowest order in as- In the frame where the vir- 
tual photon momentum is g = (z^, 0±, Vz/^ + Q'^), its light-cone components 
q^ = q^ ± q^ are 

q^ ~ 2z/, q^ ~ ~ —mx. (5) 

The Fourier transform in (|^) then implies a resolution in coordinate space of 
6y~ ~ l/2z/, ^ l/mx. In the scaling limit Q^, u —>■ 00 with a; fixed, the 
most relevant separations y between the photon currents in (^) are light-like 
distances of order the 'loffe length' l/2mx [^. 

A visualization of ImT^^^ in coordinate space is given in Fig. 1. In the 
target rest frame, where P = (m,0), the target is moving along the axis, 
while the photon enters along the positive light cone. The imaginary part 
of Tfj_u measures the overlap of two Fock states of the target, which have 
identical parton content except for the quark struck by the photon, whose 
y^ coordinate is offset by 6y^ ~ l/mx. The magnitude of the overlap is a 
measure of the mobility of the struck quark in the target wave function. 

The precise meaning of the Fourier transform of the quark structure func- 
tion is given by the operator product expansion through the relation [0, ^, Q 
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which here is formulated in A" = gauge at a renormalization scale /i^ 
and in the rest frame of the target. The left-hand side of (M) measures [cf. 



Fig. 1) the interference of target Fock states after either a quark is displaced 
along the light cone from to = y° + y^, or an antiquark is displaced 
the opposite distance. The relative minus sign is due to Fermi statistics. 
The kinematic upper limit of the x-integral on the right-hand side of @ is 
the atomic number A of the target, due to the scale m = m^v used in the 
definition of x. 

Experimental information on A-dependence is available mainly for the F2 
structure function. Subtracting from (^) the same relation with — > —y~^ 
we have 

= Aim J dx[q{x, jj."^) + q{x, jj,^)] sin ^ — '2~^ ' ^'^^ 

After summing over the quark flavors weighted by ef, the integral can be 
evaluated using the measured F2 structure function for a range of targets. 

2. Numerical Analysis 

The nuclear target effects on F2 can be expressed through the ratio 



of the structure function measured on a nuclear target A to that on the 
deuteron D. Experiments show that the ratio (^) is practically independent 
of There is data on Ra for A = Re, C and Ca from both CERN and 
SLAC [0, which together cover the ranges 0.0035 < x < 0.88, 
0.00015 < X < 0.8 and 0.0035 < x < 0.8, respectively. Measurements on heav- 
ier nuclei have been done for different choices of A in different experiments, 
and will not be used here. The available data on Rc is shown in Fig. 2. 
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In our evaluations of the Fourier transform (0) we fitted the data with a 
smooth curve (sohd hne). We also used an integration algorithm based on 
the discrete data points to obtain an error estimate. 

We used the fit of the structure function given by the NMC Col- 
laboration |ll|], evaluated at = 5 GeV^ (a typical value in experimental 
determinations of Ra/d{x)). This fit also includes SLAG and BCDMS data. 
The F^(x,Q'^ = 5 GeV^) structure functions are obtained as a product of 
the NMC fit for and our fit of the ratio The Fourier transform (|^) 
then gives the quark 'mobility' distribution in coordinate space for various 
nuclei, 

Q\y\Q') = ['dxF,^{x,Q')R^{x f''^'^''y^^^\ (9) 

JO X 

Using the quark mobility distributions Q'^iy^ ■,Q'^) we can then form the 
target ratio in coordinate space, 

Strictly speaking, the upper limit of the x-integral in should be at 
X = A. The large x region is, however, unimportant in the Fourier transform 
due to the small size of in this region. This is illustrated in Fig. 3, 
which shows the contribution to the integral in from various regions of x. 
The fact that the large x region is insignificant also implies that the effects of 
nuclear Fermi motion are suppressed in coordinate space (at moderate values 
oiy+). 

From Fig. 3 one can already anticipate important cancellations of the 
A-dependence in the Fourier integral. For y^ = y^ + y^ below 2 fm, the 'anti- 
shadowing' nuclear enhancement region beats against the EMC suppression 
region, while for y~^ around 5 fm anti-shadowing tends to offset the shadowing 
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suppression at low x. At large light-cone distances the low x region dominates 
the integral, and the shadowing effect will prevail. 

Assuming a power law dependence -F^(a;) ^ cax~^ for small x gives 



4-\ /3 
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If the power f3 is independent of A, i.e., if -R^la^) in (i) approaches the 
constant 2ca/Acd for a; — 0, the coordinate space ratio RAiy'^) in ([T0|) will 
approach the same constant 2cyi/Aci5 for — oo. This limiting value is 
numerically reached only for very large y^, however (see Fig. 4 below). 

In our numerical calculation, we have set F2{x,Q'^) = F2{x = 0.01, Q^) 
for X < 0.01, motivated by the x range of the NMC parametrization for F.^. 
We have assumed the ratio Ra{x) to approach a constant for x — 0, and 
verified that the results presented here are insensitive to the precise value 
of Ra{0) and to the behaviour of F.^ at x < 0.01, where HERA data [12 
actually show a rapid rise of F^ with decreasing x. 

3. Results and Discussion 

The coordinate space target ratio RA{y'^,Q'^ = 5 GeV^) (|10D is shown in 
Fig. 4 for y4 = He, C and Ca. Even for the heaviest (Ca) nucleus, the ratio 
is within 2 % of unity for y^ = y^ + y^ < 5 fm. At larger distances, the 
shadowing effect is clearly visible. 

Fig. 5 shows error estimates for Rc{y^,Q'^ = 5 GeV^) derived from the 
actual data points and errors of the measured momentum space ratios Rc{x) 



by using trapezoidal rule integration |jT3|. For large j/"*", the finite spacing 



of the data points in x makes any error estimate dependent on how smooth 
Rc{x) is assumed to be. The result for Rc{y^,Q'^) shown in Fig. 4 should 
thus be reliable at least for ?/+ < 7 fm. We have also verified that our results 
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are insensitive to the value of used in evaluating the structure function 
of the deuteron. 

At low , the quark mobility distribution Q^{y^) is readily seen from 
(H) to be proportional to y"*", 

Q^(y+^0)^y+-j^ dxF^ix), (12) 

where the integral measures the total fraction of target momentum carried 
by quarks. Previous careful estimates of the A-dependence of this fraction 
[0], which took into account finite energy effects in the data, gave results 
compatible with no nuclear dependence, 

dxlF^"" - F^] = (-1.5 ± 0.4 ± 1.4) ■ 10~^ (13) 

JO. 0035 

to be normalized by / dxF-P ~ 0.15. As seen from Fig. 4, the deviation 
of Rcaiu^) from unity is no larger than it is at y^ = in the whole region 
y'^ ^ 5 fm. Such a weak A-dependence of the structure function in coordinate 
space appears rather surprising, given that the EMC effect is of (9(10 . . . 15%) 
for C and Ca nuclei in momentum space. As discussed above, the weakness of 
the nuclear effect in coordinate space is due to cancellations between nuclear 
enhancement and suppression regions in momentum space. 

Some of the models discussed in the context of nuclear effects postulated 
that the effective radius of a bound nucleon would be larger by up to 30% 
compared to that of a free nucleon ("nucleon swelling" |T^). The Fourier 
transform @ provides the exact relation between distributions in coordinate 
and momentum space. We find that a simple assumption Q^^^y^) oc Q^{^y^) 
leads to an incorrect shape of the ratio RA{y^) for y^ < 3 fm, as shown in 
Fig. 6, unless 1 — ^ ^ 2%. Such a model in any case fails to describe the 
shadowing region of large y~^. 
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This analysis is not directly applicable to the "Q^ rescaling" models [|15 
where a nuclear effect arises because the effective value of is taken to be 
different for bound and free nucleons due to their different radii. 

Fig. 4 suggests that ignoring a possible ^ 2% effect for < 5 fm, the main 
nuclear effect in coordinate space is the shadowing phenomenon at large . 
The physical reason for shadowing in DIS is well understood at a qualitative 
level, and quantitative models have been successfully constructed [|1^]. It 
may thus be of some interest to see how much of the structure in momentum 
space can be ascribed solely to shadowing in coordinate space. This can be 
studied be evaluating the transform inverse to the one in (|^), 

F^[x) = — x dy^Q'^{y^,w)sm{mxy^/2) (14) 

TC Jo 

using a function Q'^{y~^,w) defined as 

^ '^^"1 y^>w, ^^^^ 

where w = (9(5 fm). Substituting the definition (^ of Q^'^^iy^) and in- 
terchanging the order of integration one obtains an expression suitable for 
numerical evaluation, 

F^ix, w) = F^{x) + - dx' iQ^ix') - g^(x')" 

TT JO L J 

{sin[m(x' + x)w/2] sin[m(x' — x)w/2] 1 
x' + X x' — X J 

The corresponding momentum-space target ratio 



(16) 



^-^^"'"^= (A/2)Ff(x) ^'^^ 

for two choices of w is compared with Ra{x) = Ra{x,w = 0) in Fig. 7. 
It can be seen that the nuclear effects observed in the data for x ^ 0.4 
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(which includes the anti-shadowing enhancement and the beginning of the 
EMC suppression) can be obtained by assuming no other nuclear effect in 
coordinate space than the shadowing for y+ > 5 fm. The nuclear effects at 
larger values of x do, however, depend sensitively also on the small (< 2%) 
effects suggested by the data at light-cone distances y'^ ^ 5 fm. 

4. Summary 

We have studied the nuclear effects on parton distributions in coordinate 
space by Fourier transforming the measured momentum space (x-Qj) distri- 
butions. Parton distributions in coordinate space can be rigorously defined 
using the Operator Product Expansion. Intuitively, they measure parton 
mobility in the target wave function, in terms of an average overlap between 
wave function components where one parton has been offset the given dis- 
tance along the hght cone. 

The parton distributions at large light-cone distances are sensitive 
only to the momentum distribution at small xbj, and thus reflect the well- 
known shadowing phenomenon. For < 5 fm {i.e., for = < 2.5 fm), on 
the other hand, we found the nuclear effects to be surprisingly weak (< 2% 
for A = He, C and Ca). Numerically, this is due to cancellations in the 
Fourier integral between regions of nuclear enhancement (anti-shadowing) 
and suppression (shadowing, EMC effect). Conversely, the observed nuclear 
dependence in momentum space reflects the mixing of effects from long and 
short distances. If the small nuclear effect at moderate hght-cone distances 
is interpreted as an effective 'swelling' of the nucleons in the nucleus, only a 
< 2% increase in the nucleon radius can be accomodated. 

Through an inverse transform we verified that the nuclear effects ob- 
served for xb] ~ 0.4 can be obtained solely from the shadowing phenomenon 
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at 7/"*" ^ 5 fm in coordinate space. The effects seen at larger xbj depend sen- 
sitively also on weak nuclear dependencies at y+ < 5 fm. 
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FIGURE CAPTIONS 



Fig. 1. A space-time picture of the imaginary part of the matrix element 
T^u of (^) for forward virtual photon-nucleon scattering. 

Fig. 2. Our fit of the structure function ratio Rc{x) for carbon compared 
with data. Squares: NMC data 0; diamonds: SLAC-E139 data [p!0|] . 
The sohd curve is the result of a least-squares cubic spline fit with 
interior knots at a; = 0.002, 0.02, 0.2 and 0.75 and an extra "data point" 
Rc{x = 0) = 0.85 added at x = 0, with error ARc{x = 0) = 0.01. The 
dashed curves show the effect of varying Rc{x = 0) between 0.7 and 
0.95. (a) Logarithmic horizontal scale, (b) linear horizontal scale. 

Fig. 3. The solid curve shows the mobility distribution Q^{y~^, = 5 GeV^) 
for deuterium. The dotted curves show the contributions to Q^{y~^, Q'^) 
from different regions in the x integral of (^. 

Fig. 4. The coordinate space ratio RA{y^, = 5 GeV^) for A = He, C and 
Ca as obtained by using cubic-spline fits of the Ra{x) data (c/. Fig. 2 
for carbon) and the NMC parametrization of F^{x,Q^) [|ll[]. At very 
large y~^, RA{y^) tends to the constant value Ra{,x = 0), which has 
here been set to 0.92, 0.85 and 0.75 for He, C and Ca, respectively. 

Fig. 5. An estimated error band for Rc{y+,Q^ = 5 GeV^). The solid line 
is as in Fig. 4. The dotted lines show error estimates based on the 



experimental errors of Rc{x) and trapezoidal rule integration |jT3|. At 
y^ > 7 fm, the trapezoidal rule integration becomes unreliable because 
of the finite spacing of the data points in x. 
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Fig. 6. Solid curve: Rciy~^,Q'^ = 5 GeV^) from Fig. 4. Dashed curves: 
Rc{y^) as obtained by setting Qc{y^) = Qoi^V^)- 

Fig. 7. The momentum space ratio Rc{x^w^Q'^ = 5 GeV^) obtained from 
the inverse transform (|T6D using the modified coordinate space distri- 
bution QA{y^,w) ([ISD, in which all nuclear effects are eliminated for 
< w. 



13 



0.0 



Mill 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


- 


Q (y ,Q =5 GeV ) 

/ ^ 

/ \ *— s 

/ \ 1 
1 ^ ■ X / 

^^-^-^ ^ ' 




^^^^ 
— / 

/ / 
/ / 
/ / 


- 


0<x<l / 

/ ' / ^ - 




/ / / \ 0.075<x<0.25- 
/ / ' ' 




rzzz": ^ ^....^ ^.^..^^1,^^ 

0.8<x<l \ / \ / 





0.25<x<0.8 V 




0.1 1 10 




[fm] 



1.1 




1.04 



1 1 — I — I — r 



1 1 1 1 1 — I — r 



Error estimate of Rc(y^) 



1.02 



1.00 



0.98 



0.96 



J I I I L 



J I I I I I L 



0.5 



1 



[fm] 



1.1 



1 — I I I I 



1 1 — I — I I I I 



1 — n 



1.0 



Rc(y^Q'=5 GeV") 



^=0.99 




0.9 



^=0.8 / 



0.8 



J I I I I 



0.5 



1 



J I I I I I I 



10 



[fm] 



20 



J u 



